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Abstract: Jet shape observables which involve measurements restricted to a part 
of phase space are sensitive to multiplication of soft gluon with large relative angles 
and give rise to specific single logarithmically enhanced (SL) terms (non-global logs). 
We consider associated distributions in two variables which combine measurement of 
a jet shape V in the whole phase space (global) and that of the transverse energy flow 
away from the jet direction, E out (non-global). We show that associated distributions 
factorize into the global distribution in V and a factor that takes into account SL 
contributions from multi-gluon "hedgehog" configurations in all orders. The latter 
is the same that describes the single-variable E out distribution, but evaluated at a 
rescaled energy VQ. 
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1. Introduction 

In a theory with dimensionless coupling one would naively expect cross sections to 
scale in a simple manner at very large Q 2 , Q 2 ^> m 2 , with m a generic mass scale. 
However, in QCD this is not true for two reasons. Firstly, due to UV divergences, the 
effective interaction strength does vary with scale (as in any quantum field theory 
with dimensionless coupling). As a result perturbative (PT) corrections to cross sec- 
tions slowly decrease with Q 2 as powers of a s oc 1/ ln(Q 2 / Aq CD ) . Secondly, there are, 
generally speaking, collinear and infrared divergences (collinear 2-parton splittings, 
soft gluon radiation). As well known (Bloch-Nordsieck [1]), in "good" inclusive ob- 
servables which do not include observation (fixing the momentum) of a single hadron, 
either in the initial or final state, the logarithmic collinear and infra-red divergences 
cancel. 

The classical examples are the total cross section of e + e~ annihilation into 
hadrons, r — > hadrons decay width. Here the collinear and infra-red divergences, 
present in both real and virtual corrections, cancel completely in the (unrestricted) 
sum. Since here we don't have dimensional parameters other than Q 2 (in the 
m 2 /Q 2 — > limit), the total cross section is given by the simple Born expression 
modulo calculable a™(Q 2 ) corrections. 



Born 



(Q) • g (a s (Q 2 )) , g(0) = 1 (Q 2 »m 2 ). 
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Moving to less inclusive measurements one faces different situations. The first 
case involves fixing (measuring) momentum of a hadron, e.g. that of the initial pro- 
ton in DIS, DY, . . . (structure functions, SF) or of a final hadron (fragmentation 
function). Then soft divergences still cancel but collinear ones do not, making such 
observables not calculable at the parton level. These effects, however, turn out 
to be universal and, given a proper technical treatment, can be factored out as 
non-perturbative (NP) inputs. What remains under control then is only the Q 2 - 
dependence (scaling violation pattern). 

Collinear divergences in the final state may be avoided altogether if one looks at 
energy flows rather than individual hadrons. More generally, one can study a family 
of the so-called collinear and infra-red safe (CIS) jet shapes V, 

V = Y J <hl (1.1) 

i 

where the sum runs over all particles (hadrons) in the final state, v(k) being a con- 
tribution of a single particle with 4-momentum k. Being (by construction) linear in 
particle momenta, such observables are also free from collinear (and soft) divergences. 

However, here the cancellation of real and virtual effects is not complete and 
leaves trace in the PT-calculable distributions over V. Indeed, taking the value of a 
generic jet shape observable 7 C 1 we squeeze the phase space thus inhibiting real 
parton production and multiplication. Since the virtual PT radiative contributions 
remain unrestricted, the divergences do cancel but produce finite but logarithmically 
enhanced leftovers: 

X(Q,V)= f dV^± = f{a s (Q 2 ),\nV). (1.2) 
Jo <?tot dv 

Each gluon emission brings in at most two logarithms (one of collinear, an- 
other of infra-red origin). These leading contributions are due to multiple soft gluon 
bremsstrahlung off the primary hard partons that form the underlying event, which 
can be looked upon as independent gluon radiation and can be easily treated. A 
clever reshuffling of PT series, based on universal nature of soft and collinear radia- 
tion (factorization) results in the exponentiated answer in the form 

oo 

In E(Q, V) = <(Q 2 ) i A n ln" +1 V + B n \n n V -\ ) . (1.3) 

n=l 

In what follows we will discuss only these two series of term and neglect subleading 

small corrections of the order of a s , a 2 \nV, The A n series is referred to as 

double logarithmic (DL) and B n as single logarithmic (SL). 

The whole leading series A n actually originates from the effect of running of 
the coupling in the basic one gluon radiation term, a s In 2 V, n — 1. Indeed, the 
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expansion in a s (Q 2 ) is pretty artificial since in reality it is a wide range of scales, 
(VQ) 2 <C k 2 <C Q 2 , at which the coupling a s {k 2 ) actually enters: 



A ias (Q 2 )\n 2 V 




: s (k 2 ) In f = £ 4» • a n s {Q 2 ) ln" +1 V, (1.4) 



* n=l 



where fct is the gluon transverse momentum and the logarithmic factor is due to the 
soft enhancement. As a result, the coefficients A n are straightforward to obtain. It 
is important to realise that what makes the coupling run in Minkowskian observ- 
ables is collinear gluon splittings (into gluons with energies of the same order, "hard 
splitting", and qq pairs) and the inclusive CIS nature of the observable (see, e.g. [2]). 

Let us remark in passing that the scales below (VQ) 2 with k 2 running down to 
zero are actually present as well. They, however, do not contribute at the PT level 
but give rise to the NP power suppressed corrections (see [2,3]). 

The subleading series B n start from the first SL correction aglnV (n = 1) (of 
either collinear or soft origin). In higher orders, n > 2, taking care of SL terms 
involves a careful treatment of a s (physical scheme), of its running argument, as well 
as precise fixing of the scales of the leading logarithmic terms in (1.3). 

In the case of the so-called "global" observables, that is measurements in which 
the observable (1.1) accumulates contributions from final state particles in the whole 
phase space, all SL contributions are generated by a careful treatment of gluon 
bremsstrahlung off the primary partons (Sudakov exponentiation). Obviously, there 
exists also gluon multiplication and, in particular, gluon bremsstrahlung off sec- 
ondary partons. As already stated, collinear "hard" gluon decays make the coupling 
run. As for soft gluon bremsstrahlung off secondary gluons, it does not affect global 
observables at the DL+SL accuracy (1.3) (the first correction being O (a 2 In V)). 
Global observables have been intensively studied in the literature both perturba- 
tively [4,5] (DL+SL as in (1.3)) and with account of the first leading NP power 
correction 0(1/ Q) [3]. 

However, as was recently noticed by Dasgupta and Salam, certain jet shape 
measurements turn out to be sensitive to soft gluon multiplication effects at the SL 
level starting from n = 2. These observables involve measurements restricted to a part 
of phase space and were correspondingly dubbed non-global [6]. The simplest example 
is given by particle transverse energy flow E ont in e + e~ annihilation, measured in an 
angular (pseudorapidity) region C out which is away from the jet direction (thrust 
axis) by a finite angle 9q: 



It is a CIS observable. Moreover, there is obviously no collinear enhancement (A n = 
0). The leading effect is SL - exponentiation of independent large angle soft gluon 




(1.5) 
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emission off the primary jets. Apart from it, however, additional contributions of 
the same order emerge here due to soft gluon-gluon correlations, which were absent 
in global observables. 

Imagine a system of gluons whose energies are strongly ordered: 

k e < < . . . < h < p ~ Q/2. 

Then only the hardest of the gluons belonging to C ont would contribute to the observ- 
able while the contributions of all other (much softer) gluons would cancel against 
corresponding virtual corrections. This means that it suffices to consider multi-gluon 
systems with the softest offspring k n e C ont being the only one to be measured, while 
all harder companions do not contribute to the observable, hi e Ci n , i < n, where 
C m is the complementary angular region, close to the jets. 

Such configurations of n > 2 gluons contribute to the integrated distribution 
(1.5) at the SL level as (a s \n(Q / E out )) n . This is the origin of the Dasgupta-Salam 
discovery. 

These specific subleading contributions are not easy to analyse analytically order 
by order. In spite of the fact that in the strong energy ordered kinematics all (n + 2)! 
amplitudes are known and given by soft insertion rules, to give a compact answer to 
the square of the n-gluon matrix element is possible only in the large- N c limit (see [7]). 
Moreover, angles between gluons are of the order one (hedgehog configurations), so 
that the phase space integrations can be handled only numerically. As we shall 
show below, the all-order resummation of these contributions can be carried out in 
the large- N c limit. This leads to an evolution equation [8] which has a highly non- 
linear structure. Therefore, its solution can be given in an analytic form only in an 
academic high-Q limit (a s \n(Q/E out ) ^> 1). 

Recall that non-linear evolution equations for generating functionals describing 
multi-parton ensembles have a long history. In the leading collinear approximation 
they, in particular, form the basis for Monte Carlo generators and for the standard 
theoretical jet studies (jet rates, in-jet particle multiplicities and spectra, etc.) [7,9]. 
When large angle soft gluon radiation becomes important (as, e.g., in inter-jet particle 
flows - the bread for the so-called string effect), in the collinear approximation (that 
is when other secondary gluons stay quasi-parallel to the primary hard parton) the 
QCD coherence is at work. As a result, a single soft gluon emission at large angles 
is very simple and is determined by the total colour charge of the jet. 

However, when the large-angle hedgehogs are considered, the new evolution equa- 
tion for such systems has an essentially different (and more complicated) structure 
which, as mentioned above, can be practically handled only in the large- N c approx- 
imation. 

It is important to stress that most of the actual experimental measurements are 
subject to non-global effects. Experiments often (if not always) involve phase space 
restriction. For example, 
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• the very first CIS jet cross section invented by Sterman and Weinberg back in 
1979 [10] provides a perfect example of a non-global observable; 

• in hadron-hadron interactions accompanying hadron production is studied in a 
limited rapidity range (in particular, the famous pedestal distributions in hard 
hadron-hadron collisions) ; 

• direct photon studies necessarily involve photon isolation criteria; 

• a family of QCD string/drag effects deals with particle flows in restricted inter- 
jet angular regions; 

• profiles of a separate jet (rather than studying shape variables of the whole 
event). For example, characteristics of the current fragmentation jet in DIS 
which is based in a one-hemisphere particle selection. 

All these and many other similar measurements contain non-global PT corrections. 
With exception of the current jet in DIS (see [11]) and the case study of the E ont dis- 
tribution in e + e~ [8,12,13], the non-global effects have not been studied theoretically 
even at level of the first (a s log) 2 correction. 

In spite of being formally subleading, these effects may significantly modify the 
PT predictions for the non-global observables, as was shown by the case study of the 
E out distribution in e + e~ by Dasgupta and Salam [12]. 

Berger, Kucs and Sterman [14] have recently formulated a programme of how 
to avoid non-global logarithms in a measurement of transverse energy flow away 
from jets (E ont ). They suggested to squeeze the jets and thus suppress multi-gluon 
hedgehogs in the C- m region. In e + e~ this amounts to introducing the associated 
distribution in two variables, E out and V (for example, V — 1 — T with T the 
thrust), and considering the region V which selects 2-jet-like configurations, 

^2ng(Q, V, E out ). 

They treated the region in which the two characteristic scales of the problem are 
comparable, VQ/E out = 0(1), so that 

a s ln — — (1.6) 

^out 

amounts to a negligible correction 0(a s ). The authors stated in the Conclusions 
to [14] that in this approximation their "formalism is sensitive mainly to radiation 
stemming directly from primary hard scattering" . 

In this paper we consider a general case 1 of E out being potentially much smaller 
than VQ. We show that, having extracted the global DL and SL enhanced term 

1 We take V to be a global observable, though one can equally well restrict V to C; n as was done 
in [14]. 
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in lnl^, the remaining SL corrections (a s \nV) n and (a s \n(Q / E out )) n conspire to 
produce the powers in (1.6). 

We analyse these non-global logarithms in all orders and demonstrate that the 
associated distributions factorize as 

WQ, V, E out) = V) ■ Knt(VQ, E out ). (1.7) 

Here H(Q,V) is the standard global distribution (1.2) and E out (Q, E out ) is the SL 
distribution in E out (at the same total energy Q) which takes into account contri- 
butions from multi-gluon hedgehog configurations in all orders and contains the full 
dependence on the geometry of the measurement (9 ). This is the main result of the 
paper. 

We will also discuss the asymptotic behaviour which according to (1.7) reduces 
to that of the simplest non-global distribution S out which had been studied in [8]. 

Let us stress that one needs to keep under the best possible control effects in- 
duced by radiation of relatively soft gluons not merely for the sake of improving PT 
predictions. What makes such studies even more important and interesting is the 
fact that the physics of small transverse momentum gluons is that of confinement. 2 



2. Observables 

In this paper we will consider e + e~ annihilation into hadrons and define a cone 
around the thrust axis and denote by C out and C in the regions outside and inside the 
solid cone (inter and intra jet regions). 

We start by considering the first factor in (1.7) that is the (integrated) global 
distribution H{Q,V) defined as 

E(Q, V)=J2 / ^© (v- £ v(k t )\ , (2.1) 

with the sum extended to the full phase space C tot = C in + C out . Here da n and a tot 
are the production and total cross sections. The probing function v(k), linear in the 
particle momentum k, can be represented as 

v(k) = | • h( V ) , (2.2) 

with k t the transverse momentum with respect to the thrust axis and i] the pseu- 
dorapidity. To get the distribution in thrust, T, we have to take h(rf) = and 
the set V = 1 — T. For small 1 — T, summing over final particles we may neglect 
contributions from the primary hard partons (quarks) since h(r)) vanishes at large 
1 77 1, and restrict ourselves to to considering only secondary soft gluons. 

2 In market terms, a $10 6 problem [15] 
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The distribution in broadening B, for example, corresponds to h{rj) = 1, V = 2B. 
In this case, even for B <C 1, one needs to include the contributions from the recoiling 
primary quarks. 

Now we define analogously the non-global distribution (the second factor on the 
r.h.s. of (1.7)): 

£ ou t(Q,£out) = J2f ^ ~ k A ' ( 2 - 3a ) 

where we have chosen to measure the transverse energy flow accumulated in the 
angular region C out away from the jets. 

Finally, as was suggested by Berger, Kucs and Sterman [14], we introduce the 
shape observable distribution in two variables (the associated distribution) 

z 2ng (Q,v,E out ) = Yl I — ®(v- Yl ^fcOJ-efsout- k A- ( 2 - 3b ) 

n J °" tot V ieCtot / V iec out J 

3. Analysis and resummation 

We are interested in the soft region 

V,£ ou t«Q. (3.1) 

3.1 Mellin factorization 

First we recall how the resummation with the SL accuracy is done for global ob- 
servables. To this end one has to factorize the theta- function in (2.1) by the Mellin 
transform: 



£«,V)=/^e-E W ,,), 

t(Q, V )=y, j — n "<*.). "<*)= e 



da ^ ( 3 - 2 ) 



As has been mentioned above, in the thrust case the contributions of the primary 
partons can be dropped (since v(k) vanishes at large \rj\). 

In the case of broadening the recoiling primary partons do contribute. This 
causes certain complication. This recoil effects can be taken care of by expressing 
the momenta of recoiling quarks in terms of those of the secondary partons via 3- 
momentum conservation. Thus, also in the case of broadening one can restrict the 
product in (3.2) to the secondary soft partons, provided one factorizes also the mo- 
mentum conservation constraints by the Fourier transform. This amounts to incor- 
porating into £ an additional Fourier variable associated with transverse momentum 
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conservation. Peculiarities of B as a variable employed to squeeze the jets, can be 
analyzed along the lines of [5]. This leads to the known modification of the global 
distribution but will not affect the dependence on E ont in (1.7). 

In what follows to simplify the discussion we neglect this complication and con- 
centrate on the thrust case. 

For the out-distribution we have 

£out(Q,£out)= / 4^e"° UtEoUt/Q Zout(Q,Vout), 

J Zinv out 

S ou t(Q,^out) = ^ / " H M °ut(^), ( 3 - 3 ) 

m J °~tot -rr< 
n «GGtot 

u out (k) = 0- m + 9 out e-^ k */Q = 1 - 9 out (1 - e -"»tfc«/Q) . 
Here 9 in , 9 out are the support functions for k in the regions C m and C ou t, respectively: 

in = tf(M-77o), ^out = #(vo-\v\), (3-4) 

with 770 the pseudorapidity corresponding to the opening (half-) angle of the cone, 

i] = -In tan— . 

Note that the fact the probing function u out (k) — 1 for k G C in means that this 
parton does not contribute to the measurement. 

Finally, for the associated distribution we introduce a double Mellin transform: 



WQ, V, E ont ) = [ dU dVw * e ^t WO s 2ng (Q, u, u out ) , 
J Zixi v Zixi u out 

S 2 ng(Q, V, Vaut) = ^2 " J I U 2ng(fci) ■ 



(3.5) 

''2ng(fcj) • 

ieCtot 

The new source function is given by the expression 
W2ng(£0 = e-" v W (9 in + ^oute _I/outfet/Q ) = e- vv ^-6 mst e- vv ^ (l - e"" -***/ ) . (3.6) 

It corresponds to measuring V everywhere and E ont only in C out . 

In the Mellin space all the distributions become exponents of the so-called radi- 
ators: 

t(Q,v) = e- R «>"), (3.7) 
t oat (Q,u out ) = e- R ^ {Q ^\ (3.8) 
S 2ng (g, !/, z/ out ) = e --R**«."."<»t) . (3.9) 
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3.2 Global radiator 

For the global distribution the answer is simple and the radiator is known to be given 
by the SL-improved Sudakov expression 

R(Q,v) = 2C F [ Q ^^l f^dr, (Is-™ <*>), W = ln^. (3.10) 
Jo 71 Jo 

In the definition of rj milx we have incorporated a SL correction to the "hard" piece of 
the quark splitting function 

l + (l-a) 2 ( a\ k t _ 

' a = — e '. 



= g = I 2 J ' 



It is straightforward to verify then that the virtual SL contribution due to large ra- 
pidities, a ~ 1, can be embodied into rescaling of the upper limit of the rj integration 
as follows: 



/ daP{a) =^ [ 

Jkt/Q JO 



^?7max 

dr\. 

'k t /Q 



The two-loop radiator (3.10) (with the coupling a s in the physical "bremsstrahlung" 
scheme [16]) has single-logarithmic accuracy. This means that the neglected correc- 
tion is O (ag In i/) which translates into \og(Q/V) (see, e.g. 



3.3 Non-global radiator 

The radiators for both out and associated non-global distributions, (3.8) and (3.9), 
have the structure 

R A = rW + R%\ A = out, 2ng . (3.11) 

Here R^ is the standard Sudakov one-gluon emission contribution and is the 
term due to correlated multi-gluon radiation. (The latter is absent (negligible within 
the SL accuracy) in the case of a global observable.) 

Sudakov piece. The Sudakov piece R A \ 

i#> = 2C F [ ^^Ml P ax dr)[1 _ UA{k) ] , (3.12) 

J Kki- TT Jo 

is given by the expression (3.10) with proper probing functions: 

[ 1 - UoutW ] = ^out (1 - e-^*/ Q ) , (3.13a) 
[ 1 - u 2ng (k) ] = [ 1 - u{k) ] + e out e~^ (1 - e-^/Q) 

= [ 1 - u(k) ] + [ 1 - Uout (k) ] . (3.13b) 

We compute first R^t- Substituting (3.13a) into (3.12) we obtain 

lC(Q, "cut) = 2C F [ Q ^^Ml P dr] [i_ e -*»*M>] , (3.14) 
Jo K t 71 Jo 
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with the rapidity integral restricted to the C out region. Here we have used that 
i]o ~ 1 < ?7 m ax- This expression does not have a collinear singularity and therefore is 
a SL function. Therefore we can approximate [1 — u out ] in (3.14) by a theta-function 
$(k t — Q/v out ) and evaluate the inverse Mellin integral by simply substituting Q/E ont 
for u ont . We arrive at 

tfout(^out) ^77o-A(g,Kut) = 1lV(Q,E out ), (3.15) 

where we have introduced a convenient SL function 

AiQ , E)=C J Qd -tl^l. (3.16) 

J E "'t 7T 

Similarly, using (3.13b) in (3.12) gives the expression for the Sudakov piece of the 
associated distribution radiator: 

R% g (Q,is,is ut) = R {l \Q,v) + 5R {1) (Q,v,v out ). (3.17) 

Here the first term reconstructs the global distribution Y,(Q,V) in (1.7) and the 
addition contribution 5R^ to the Sudakov radiator reads 

5R% g (Q, u, VwA ) = 2C F [ Q ^^Ml P dr] e — <*) (l _ e -w (318) 

J ti t TT Jq 

This expression is identical to that for the out-case, (3.14), apart from the additional 
factor exp(— v v(k)). Since (3.18) is a subleading SL correction, we can treat this 
factor, within our accuracy, as simply proving an additional restriction upon the 
transverse momentum. As long as i] = O (1), we have v{k) ~ k t and thus 

e -v(k)u ^ &(Q-vkt). 
In V^-space this translates into the condition 

k t < VQ < Q. 

Then the integral for 5R^ becomes the same as that for R^ t with Q replaced by VQ: 
5R^ g (Q,u,u out ) =*► ^ m -A(VQ,E out ) = TZ^ t (VQ, E out ). (3.19) 



We conclude at the level of the Sudakov contributions to the radiators the factorized 
answer (1.7) holds. 

Now we have to analyse the specific correlation contribution R^ in (3.11). This 
is also a SL function whose PT expansion starts at the level of (a s log) 2 . 
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2-Loop correlation. At two loops is given by the integral 



= J du 2 (h, k 2 ) [ U A (h, k 2 ) - u A (k 1 )u A (k 2 ) } , (3.20) 

where du 2 is the distribution for the correlated emission of two secondary partons off 
the primary qq system. Here U(k) is an "inclusive" source function for the parent 
massive gluon. As was shown in [17], to properly reconstruct the two-loop Sudakov 
expression with the running physical coupling, this source has to be defined in terms 
of the single parton source function UA(k) = UA(k t ,rf) as 

U A {ki,k 2 ) = u A (Vkt +m 2 ,r))] k = h + k 2 , m 2 = k 2 . (3.21) 

Here rj is the rapidity of the parent, 



V^ 2 + m 2 e ±r > = k tl e ±rn + k t2 e ±m . 

The distribution duo 2 is singular at m 2 = when the two partons become collinear 
or one of them is much softer than the other. In particular, in the collinear limit 
the secondary partons obviously belong to the same angular region {C m or C ut) 
and then, as can be easily seen from (3.20), the combination of the sources in the 
square brackets vanishes thus regularizing the singularity. The same is true for the 
soft (energy ordered) two-gluon configuration, provided the gluons are in the same 
angular region. This results in a negligible subleading contribution O (a 2 log). 

As we will see shortly, the only relevant logarithmically enhanced contribution 
emerges in the case when the two secondary partons belong to the complimentary 
angular regions, C in and C out , so that the sources u(k\) and u{k 2 ) become different 
and the cancellation gets broken. This is the configuration that gives rise to the 
non-global SL corrections, as was found by Dasgupta and Salam [6]. 

Given that there are no collinear logarithms, we can simplify (3.20) by taking 
strongly ordered parton energies, 

k t2 < hi < Q ■ (3.22) 
In this configuration the correlated 2-gluon distribution doj 2 in (3.20) reads [18] 

duj 2 (k ± , k 2 ) = 2C F CA(^yf[(^^Yc(n 1 ,n 2 ), (3.23) 

i=l ^ ' 

The function C depends on parton angles and is given by the expression 

„. . (nn) (nn) (nn) (nn) 

C(ni,n 2 ) =- — — + 



nni)(nin 2 ){n 2 n) (nn 2 )(n 2 rii)(nin) (nn x )(nin) (nn 2 )(n 2 n) (3.24) 
= WnnM \ w nni (n 2 ) + w nin (n 2 ) - w nn (n 2 ) } , 
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where 

w ab(c) = t w , v , (n a n h ) = 1 - cos a6 
{ac){cb) 

and n, n stand for the direction 4- vectors of the primary partons and n», i = 1,2 
of the secondary gluons. The angular distribution (3.24) becomes singular when the 

— * — * 

gluon momenta h\ and h 2 are parallel. 

First we observe that in the soft limit (3.20) we can put m 2 = in the source U 
to approximate 

[ U A (h, h 2 ) - UAih^UAih) ] ~ u A (h) [ 1 - u A {h 2 ) ] . (3.25) 

Let us first analyse the more complicated case of the associated distribution. Using 
the explicit expression for the source (3.6) we obtain two potential contributions: 

h\ € Qn, ^2 £ Cout : . 
n n r / 1 \ , -> (3.26a) 

e -^(fci) j-j _ e -^(fe 2 )-z,outfe 2t j ^ min{£ out , FQ} < fc 2t < fctx < VQ; 

fci £ Cout, fc 2 G C in : 
e -^(fei)-.outfeu [ X _ e -^(ta)] =^ < A; 2 , < fc tl < min{£ out , VQ}. 

We see that the only contributing region is that in (3.26a) where a harder gluon is 
close to the jet axis while a softer one contributes to the out-of-jet fc^-flow, provided 
-E'out ^ VQ. Therefore in (3.20) we can substitute 

[f/ 2 „ g (A;i,A; 2 )- M2ng (A; 1 ) M2ng (A; 2 )] 9- m (hi)0 out (k 2 ) ${VQ - h lt ) $(h 2t - E out ) (3.27) 

and we get at two loops 

R2I =► % ^ 2 (VQ, E out ) ■ F{no). (3.28) 



Here F(r) ) is given by the angular integral 

J2n r J^Cl ^2 



^Wc^/c ^(n,*)^-!*^)- (3.29) 

Oi n «/ Oout 

The same conclusion holds for the out-distribution: 

[Uoutih, k 2 ) - «out(A;i)«out(A;2)] M^i^out^) - feu) ^(Afc - E out ). (3.30) 

Comparing this with (3.27) we see that this case reduces to setting V = 1 in (3.26), 
and we obtain 

R^t =► |a 2 (Q,^)'%). (3.31) 

The result (3.28) is just the same expression that we derived for the away from jet 
distribution, except that what was Q there is now replaced by VQ. 
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This completes the proof of (1.7) with account of the 2-loop gluon correlations. 
So we can write the integrated non-global distributions (2.3) with the SL accu- 
racy as 

£ ut(Q, E out ) = e -*&«.*»t) e -^«« ) (3.32a) 
E 2ng (Q, V, E out ) = e - n o^ VQ ' E °^ e -<c«.^^t) . e( q > y )? (3. 32b ) 

where was defined in (3.15). We have shown in this section that with account 
of the 2-gluon correlation (O (A 2 )) 

K% g (Q,V,E ont ) = K^(VQ,E out ). (3.33) 

In the following section we show that (3.33) hold at all orders in A. 

4. Multi-gluon correlations in all orders 

Since there is no collinear singularities in R(°\ all SL contributions arise from en- 
sembles of soft gluons with strongly ordered energies - the leading logarithmic soft 
approximation. As a result of the strong energy ordering, it is only the hardest 
of these soft gluons in the angular region C ont that contributes to E out (since the 
contribution of softer gluons are negligible in the leading soft approximation). 

Therefore, to collect these specific SL contributions to the non-global distribu- 
tions (both the out- and associated distributions) it suffices to analyse multi-gluon 
systems with many gluons in C in and a single one in C out . Let us denote the mo- 
mentum of this gluon by q. Then, among the gluons h L G C in it suffices to consider 
only those that are harder then q, since only such gluons may affect the radiation 
of q. Softer ones, ki <C q, don't contribute to the measurement due to real-virtual 
cancellations (except as power-suppressed corrections). 

With account of multi-gluon effects the non-global correlation function be- 
comes much more involved. At the same time, the structure of its dependence on Q 
and V remains the same as in the 2-loop case. Indeed, the logarithmic integration 
in the energy of the hardest gluon k\ in the C- m region is limited from above either 
by a pure phase space restriction k w < Q for the case of TZ^t or > alternatively, by 
k\o < VQ in the case of the associated distribution TZ^ g . Since this upper limit is 
the only place where the dependence on Q (QV) enters, we automatically derive the 
relation (3.33). Therefore it suffices to consider TZ^ t (Q, E ont ). Hereafter we suppress 
the subscript and study 

rt c \Q,E ont ) = K { :l{Q,E mt ) = n%JQ,V = l,E out ). 
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Generalisation. For the sake of generalisation the 2-loop result can be rephrased 
as follows. We take k G C in and q G C out , q <C k, and write (a s = Caols/tt) 



<n(c) (n TP \ f Q dk ° f dQk - ( \ 

K (> (Q,E out ) = — ——a s w nn {n k ) 
JE out h J Cin Air 



(4.1) 

X { r 'nn k (ko, E out ) + r nkn (ko, E out ) — r nn (ko, -^out)} j 

where r ab is the 1-loop radiator for the out-distribution generated by a colorless dipole 
ab with a,b G C in : 

r a& (Q,^out) = -t^ / — / -r^ a s ^ a6 K) . (4.2) 

In the large- iV c limit the structure of soft multi-gluon multiplication can be described 
in terms of the iterative procedure (see [7]): 

W ab (l, ...,m)= w ab (£) W a£ (l, ...,£- l)W eb (£+l, . . . , m), W ab (l) = w ab (l), (4.3) 

where £ is one of m soft gluons emitted by the colour singlet dipole (ab) . We choose 
£ to be the hardest of the secondary gluons. Following this path we arrive at the 
following generalisation: 

Q 

Z a n k (ko, E out )Z nkb (ko, -Eout) 



^i?(Q>^out) = ®sW ab (n k ) jl - 



Z a b(ko, E ont ) 



(4.4a) 

where 

- In Z ab (E, E out ) = r ab (E, E out ) + TZ^ (E, E out ). (4.4b) 

In (4.4) Z ab is the generating function describing the gluon cascade which originates 
from the afr-dipole. We remark that the coupled equations (4.4) describing the large- 
angle (hedgehog) gluon configurations have a highly non-linear structure. 

Once the equation is solved, we substitute the directions of the primary partons 
nn for ab to obtain 

n { :i(Q,E out ) = TZ^(Q,E out ). (4.5) 

Recall that to obtain the all-order correlation function TZ^ g (Q, V, E out ) for the asso- 
ciated distribution we simply have to substitute Q — > VQ in (4.5). 

The asymptotic regime. The behaviour of the answer in the asymptotic limit 
a s \n(Q/E out ) ^> 1 (a s \n(VQ / E out ) ^> 1) was analysed in [8]. For the partons a and 
b inside the same C in cone and 9 ab <C 1 

Z ab ~ h , 9 C ~ A(77o)e-i A , c ~ 2.5 (4.6) 



14 



with c a universal number and A a finite factor which depends on the opening angle. 
Here 9 C is a "critical angle" which depends on the SL parameter A defined in (3.16), 
and h(x) is a universal function which is normalized, h(0) = 1, and decreases fast 
for large x as 

h(x) oc exp(— - In 2 x). 
When A is large, Z can be approximated by a theta-function 

Z ab (V,E out )~#(0 c (A)-9 ab ). 
The dependence on 9 becomes subleading in A, and we get an asymptotic estimate 

n^(V,E ont ) = c-A 2 (1 + (A- 1 )). (4.7) 

Let us remark that the TZ^t (Sudakov) contribution to InS is in this limit compa- 
rable with the subleading term in (4.7) which is geometry-dependent and could be 
computed by a numerical solution of the equation (4.4) or by the methods of [6,12]. 

5. Conclusions 

We considered the associated measurement of a global event shape variable V (e.g. 
thrust) and the away from jets transverse energy flow E out (the "flow/shape" cor- 
relation [14]) in e + e~ annihilation. We showed that the corresponding integrated 
distribution S 2ng defined in (2.3b) factorizes into the product of the global distri- 
bution in V (2.1) and the (non-global) SL distribution in E ont (2.3a) taken at the 
rescaled energy VQ, as stated in (1.7). This result holds with DL+SL accuracy, 
including resummation of the SL terms in the ratio of the two characteristic scales, 
(aMVQ/E out )) n . 

This result hold also for other shape variables, in particular, for the broadening 
B, in which case the analysis of the global spectrum is more involved [5]. 

In the most interesting limit VQ ^> E out the distribution S 2n9 describes also the 
case when the V observable is measured only inside C in (rather than in the whole 
phase space), as has been suggested in [14]. 
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